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Binary symmetry constraints are applied to the nonlinearization of spectral prob- 
lems and adjoint spectral problems into so-called binary constrained flows, which 
provide candidates for finite-dimensional Liouville integrablc Hamiltonian systems. 
The resulting constraints on the potentials of spectral problems give rise to a kind 
of involutive solutions to zero curvature equations, and thus the integrability by 
quadratures can be shown for zero curvature equations once the corresponding bi- 
nary constrained flows are proved to be integrable. The whole process to carry out 
binary symmetry constraints is called binary nonlinearization. The principal task 
of binary nonlinearization is to expose the Liouville integrability for the resulting 
binary constrained flows, which can usually be achieved as a consequence of the 
existence of hereditary recursion operators. The theory of binary nonlinearization 
is applied to the multi-wave interaction equations associated with a 4 X 4 ma- 
trix spectral problem as an illustrative example. The Backlund transformations 
resulted from symmetry constraints are given for the multi-wave interaction equa- 
tions, and thus a kind of involutive solutions is presented and the integrability by 
quadratures is shown for the multi-wave interaction equations. 



1 Introduction 

The nonlinearization technique arose in the theory of soliton equations ten 
years agouS. One gradually realizes that it provides a powerful approach for 
analyzing soliton equations, in both continuous and discrete cases, especially 
for showing the integrability by quadratures for soliton equations. The ma- 
nipulation of nonlinearization leads to finite-dimensional Liouville integrable 
systems (for example, see u~El) and also establish a bridge between infinite- 
dimer^kmal soliton equations and finite-dimensional Liouville integrable sys- 
tems ErEi WhatJSj-niore, it paves a method of separation of variables for 
soliton equations E3~t3, and thus the determination of finite-gap solutions to 
soliton equations is transformed into solving the corresponding Jacobi inver- 
sions problems PJ'Ej. 

Mathematically speaking, much excitement in the theory of nonlineariza- 
tion comes from a kind of specific symmetry constraints, engendered froni 
the variational derivative of the spectral parameter of spectral problems EJiIj, 
.which has a close relation with eigenfunctions of hereditary recursion operators 
I1ZH13. Symmetry constraints are themselves a common conceptional umbrella 
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under which .one can manipulate both mono-nonlinearizationEl and binary non- 
linearizationlij, and the study of symmetry constraints is a crucial part of the 
kernel of the mathematical theory of nonlinearization. 

Although it is not clear how to determine pairs of canonical variables to 
obtain Hamiltonian structures of constrained flows in a general case of mono- 
nonlinearization, there exists a natural way for determining symplectic struc- 
tures to exhibit Hamiltpnian forms for binary constrained flows while doing 
binary nonlinearization cB. Therefore, one can expect to establish a rigorous 
theory for binary nonlinearization. In this paper, we would only like to discuss 
the procedure of manipulating binary nonlinearization. More considerations 
on the topics related to binary nonlinearization will be made elsewhere. 

The paper is structured as follows. In section 2, a procedure for doing 
binary nonlinearization is depicted for general matrix spectral problems. Then 
in section 3, an illustrative example is carried out to make an application of 
the procedure to the multi-wave interaction equations associated with a 4 x 4 
matrix spectral problem. Finally in section a summary is given along with 
some concluding remarks. 



2 Binary Nonlinearization 

In this section, we would like to describe the procedure of manipulating binary 
nonlinearization for soliton equations (see, for example, EBLj). Let us start 
from a general matrix spectral problem 

(/> x = U(l>=U{u,X)(t>, U=(U ij )rxr, 0= (01, •••,<M T (2-1) 

where A is a spectral parameter and u = (iti, • ■ • , u q ) T is a q-dimensional vector 
potential. Associated with the spectral problem (2.1), suppose that for each 
integer m > 0, we have the evolution law for the eigenfunction <\>: 

<kn =V i - m U = V {m) {u,u x ,--- ] \)<p, FW = (^') m , (2.2) 

If the Gateaux derivative U' of U with the potential u is injective, then the 
isospectral (A* m = 0) zero curvature equation 

Ut m - V x {m) + [U, V {m) ] = (2.3) 

will determine an evolution equation 

SH 

ut m = X m (u) = JG m = J^^-, u = u(x,t m ), (2.4) 

ou 
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where X m depends on the potential u and its spatial derivatives of u up to 
some finite order, J(u) is a Hamiltonian operator, and H m (u) is a Hamilto- 
nian functional. Evidently, the compatability condition of the adjoint spectral 
problem and the adjoint associated spectral problem 

iP x = -U T (u,\)ij, ^ tm = -V^ T (u,X)^, ^ = (^---,W) T (2.5) 

is still the same as the zero curvature equation (|2.3|), guaranteed by the evo- 



lution equation u tm = X m (n) defined by (2.4). 
It has been pointed outEiJO that 

S\ , T WM f T dU(u, A) 



is a symmetry of the evolution equation ( |2.4| ), where E is called the normalized 
constant, and Q = (0, T) if u is assumed to be periodic with period T and 
51 = (-co, oo) if u is assumed to belong to the Schwartz space. Let us now 
introduce N distinct eigenvalues Ai, • • • , Ajv, and so we have 

4> { x s) = U(u,X s )^ s \ ^i s) = -U T (u,X s )^ s \ l<s<N; (2.6) 

4 S 2 = V^(u,\ s )^ s \ 4^ = -V^ T (u,\ s )^ s \ 1 < a < N; (2.7) 

where the corresponding eigenfunctions and adjoint eigenfunctions are denoted 
by (j)^ and ip( s \ 1 < s < N. Suppose that the covariant G mo is Lie point, 
and then the so-called binary Bargmann symmetry constraint reads as 

X mo = W Ms J^i, i.e., JG mo = /f^' ^^^ "! (2.8) 

s— 1 s— 1 

where 1 < s < N, are arbitrary non-zero constants, and E s , 1 < s < N, 
are iV nor mali zed constants. The right-hand side of the binary symmetry 
constraint ( [2.8j ) is a linear combination of N symmetries E s JSX s /5u, 1 < s < 
N, and the requirement of /i s ^ 0, 1 < s < iV, is natural , sin ce the symmetry 
J5X s /Su is just not involved in our symmetry constraint ( |2.8|) if some constant 
/i s = 0. It is also worthy to mention that <^ s ' and ip( s '> can not be expressed 
in terms of x, u and spatial derivatives of u to some finite order, and thus such 
symmetries (or the corresponding covariants 5X S /Su, 1 < s < N) are not Lie 
point, contact or Lie-Backlund symmetries. 



Let us now assume that we can solve the symmetry constraint (2.8) for it: 



^W,..,^);^,..,^)). (2.9) 
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Replacing u with u in the system ( fOi| ) and the system (2.7), we obtain the 
so-called spatial binary constrained flow: 

4 s ) = U(u,\ s )<t> {s \ </£ s) = ~U T (u,\ s )^ s \ 1 < s < N, (2.10) 

and the so-called temporal binary constrained flow 

= V {m) («, A s )0 (s) , Vtt 5 = -^ (m)T (w, A s )V- (s) , 1 < a < AT. (2.11) 

Note that the spatial and temporal binary constrained flows ( [2.10 ) and ( [2.11 ) 
are nonlinear with respect to eigenhmctions and adjoint eigenfunctions, but 
the original spectral problems (2.6) and (2.7) are linear. Moreover, the spatial 
binary constrained flow ( 2.10| ) is a system of ordinary differential equations, 
and the temporal binary constrained flow (2.11) is usually a system of par- 
tial differential equations but it can be transformed into a system of ordinary 
differential equations under the control of ( [2.1C ) . 

The principal task in the theory of binary nonlinearization is to show that 



the spatial binary constrained flow (2.10) and the temporal binary constrained 



flow (2.11) under the control of (2.10) are Liouville integrable, and thus the 



integrability by quadratures can be shown for evolution equations possessing 
zero curvature representations. The main tools to do this are Lax represen- 
tations generated from the stationary zero curvature equations and r-matrix 
formulations of Lax operators Ell. These will form the important part of the 
theory of binary nonlinearization and will be shed right on in our example. 
Now if </>( s ) and ij)^ s \ 1 < s < N, solve two binary constrained flows 



( 2.10 ) and ( 2.11 ) si mul taneously, then u = u will pr esen t a solution to the 
evolution equation (2.4), since the evolution equation (2.4) with the potential 
u = u is the compatability of two binary constrained flows ( 2.1 d| ) and (2.11). 
It also follows that the evolution equation u tm = X rn (u) is decomposed into 
two finite-dimensional Liouville integrable systems, and the constraint u = u 
gives rise to a Backlund transformation between the evolution equation and 
the resulting finite-dimensional Liouville integrable systems E3. Furthermore, 
two binary constrained flows (2.1C) and (2.11) can be solved by separation 
of varaibles and all what we need to do -upon obtaining separated variables 
is to solve the Jacobi inversion problems 113. This whole .process to carry out 
symmetry constraints is called binary nonlinearizationli3'L2l. 



3 Applications to Soliton Equations 

By a soliton hierarchy we means a hierarchy of evolution equations with a 
recursion relation as follows 



X m (u) = $ m X = JG m = J 



5H m 
Su 



m > 0, 



(3.1) 
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where u is the vector potential u = u(x,t m ), X m , m > 0, depend on the 
potentional u and its spatial derivatives of u up to some finite order, $ is a 
herediatry common recursion operator, J(u) is a Hamiltonian operator, and 
H m (u), m > 0, are the Hamiltonian functionals. Moreover, we have 

[X m , X n ] := (x m (u + eX n - X n {u + eX m j) = 0, m,n> 0, (3.2) 



e=0 



SH m T T 5H n 



{H m , H n }j := J (-^Y J-gf- dx = 0, m, n > 0, (3.3) 

where [•, •] and {•, -}j are called the commutator of vector fields and the Poisson 
bracket assocaited with the Hamiltonian operator J, resprectively. It follows 
that each equation ut m = X m in the soliton hierarhy ( |3.l[ ) has infinitely many 
symmetries {X n }^ =0 and infinitely many conserved functionals {-Hn}^o- This 
property can be obtained from a bi-Hamiltonian formaultion, and its natrural 
bi-Hamiltonian formaulation for the hierarchy ( |3.1[ ) is formed by J and M = 
Jfy if they are a Hamiltonian pair, where ^ denotes the adjoint operator of 
i.e., ^> = <&t. The above property can aslo be obtained from_a weaker condition 
that M is a symmetric operator (M^ = M), i.e, <&J = J^ES. 

A soliton hierarchy can usually be constructed from a hierarchy of zero 
curvature equ ations as the compatability conditions of a g iy^r -mntrW spectral 
problem {2A) and its associated spectral porblems (2^) ot3. The heredi- 



tary recursion operator can be generated from the matrix spectral problem 
), and the Hamilotnian formulations can be determined by a trace identity 
All this information shows that we can have everything for manipu- 



;.2(i 



lating binary nonlinearization for our soliton hierarchy (3.1). Applications of 
binary nonlinearization have been made to the KdV, the MKdV, the AKNiL 
e-JCaup-Newell, the Dirac, the WKI soliton hierarchies and so onEj'EErEj 
"E3. In what follows, an illustrative example of the multi-wave interaction 
equations associate with a 4 x 4 matrix spectral porblem will be carried out 
for making applications of the procedure of binary nonlinearization. 



3.1 Multi-wave Interaction Equations 

Let us start from a, Ax A matrix spectral problem: 



4> x = u (U, \)4>,U: 



/ «iA Ui2 U13 Ui4 \ 
U 2 1 a 2 A U23 «24 

U31 U32 a 3 \ m 34 
\ u 4 i U42 U43 a 4 A J 



= XU + U l7 <j> = 



h 



(3.4) 
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where Uq = diag(ai, a%, «3, 0:4), and ai, «2, 0:3, a^ are distinct constants, and 
the potential u is defined by 

U = p(Ui) = (U21,U12,W31,W13,U41,M 3 2,U23,U14,U42,U24,M43,U34) T - (3.5) 

The associated spectral problem is chosen as 



h = V (u, A) 0, y 



^ /3lA U12 U13 U14 \ 
«2l /? 2 A «23 V 2 4 
V 3 1 V32 faX W34 

\v 41 v 42 W43 AiA/ 



= XV + V 1 , 



(3.6) 



where Vo = diag(/3i, fa, 03, fa), and fa, fa, fa, fa are distinct constants. Then 
the isospectral (A* = 0) compatibility condition t/t — 14 + ft/, V] = of the 
spectral problem (|3.4[) and the associated spectral problem (3.6) gives rise to 

U u - V lx + [U u Vi] = 0, [U a , V{\ = [Vo, U{\. (3.7) 

These two matrix equations lead to the following multi-wave interaction equa- 
tions 



fa ~ fa 
iHj,t = u 



4 

. 1 rk — fa Pk ~ fa \ 1 ^ ■ 1 ■ ^ . to o\ 

i 3 .x+ } ( -)UikU k j, 1 < 1 + J < 4, (3.8) 

^ ak - ai ak - aj 

k^i,j 



which contain a couple of physically important nonlinear models as special 
reductions Ej. Note that the compatability condition of the adjoint spectral 
problem and adjoint associated spectral problem: 

i) x = -U T {u,X)i), ^t = -V T (u,\)*P, ^ = (^,^2,^8,^4) T (3.9) 

still gives rise to the above multi-wave interaction equations. It is easy to find 
that the multi-wave interaction equations (3.8) have the Hamiltonian structure 



5H 

Sua 



1 < i ^ 3 < 4, 



(3.10) 



where the Hamiltonian operator J reads as 

J = diag((ai - 0-2)00, (ai - a 3 )a 0l (ai - 0:4)00, 

(a 2 - "3)0-0, (a 2 - 04)0-0, (a 3 - 0-4)00), o = ^ 

6 



1 

-1 



(3.11) 



The Hamiltonian function H in ( |3.10 ) can directly be computed but it is 
omitted here due to its complicated form. We remark that the soliton hierarchy 
containing higher-order multi-wave interaction equations can also be generated 
from the spectral problem (3.6), but in this report we focus on the multi-wave 
interaction equations ( ft.§| ). 

3.2 Binary Constrained Flows 

Let us now construct binary constrained flows of the multi-wave interaction 
equations (3.8). Based on the general procedure of binary nonlinearization, 
the symmetry of (3.8) generated from the conserved functional X(u) reads as 

E /2l = Ep([U ,p-H 6 ^)}) = pfltTo.p- 1 ^^)]), (3-12) 



where E is the normalized constant and p is the mapping defined by (3.5). 
Upon introducing N distinct eigenvalues Ai, • ■ • , Xn, we obtain N replicas of 
Lax systems: 

4> { x s) = U(u, Xs)^ {s) , V4 S) = -U T (u, X s )ij/ s) , 1 < s < N; (3.13) 

4 S) = V(u, X s )^ s \ 4 S) = -V T (u, X s )^ (s \ l<s<N; (3.14) 

where the eigenf unctions (j)( s > and the adjoint eigenfunctions ip( s > are assumed 
to be denoted by 

= W>l S ,02s>3 S ,04s) T , ^ {S) = (As,i>2s,i>3s,lp4 S ) T , 1 < S < N. (3.15) 

These N eigenvalues lead to a more general symmetry of fl3.8p 

Z := p([U , p~\J2 ^ (s)T ^|^V)]) = pWo, £ Ps^^ s)T ]), 

s— 1 s— 1 

where the /i s 's are arbitrary non-zero constants. On the other hand, we can 
find by inspection that the multi-wave interaction equations ( |3.8| ) have a Lie 
point symmetry 

Y := p([T, E7i]), r = diag( 7l , 72 , 73 ,74), 1 < i ? j < 4, 

where 71,72,73,74 are distinct constants. In fact, we can directly prove that 
both Yq and Zq are two symmetries of the multi-wave interaction equations 



(|3_8|). Define 

N 

SU X = [T.tfi] or [U ,J2^ {S) ¥ S)T ], 



then [£/o,5Vi] = [Vb,<5?/i] uniquely determines 



N 



SVr = [F,Vi] or [F , J> s (s ¥ s)T ]- 



s=l 



Now it follows from (3/7) that the symmetry problem only requires to verify 
that 



JY 



N 



(dUuSVt) = ([r^ur,^]) or ([tf ,$>^ (s ¥ a)T ], [v ,J2^ {s) ^ s)T ]) 



satisfies the linearized system of the multi-wave interaction equations ( |3.8| ): 

(fl7i)t - (SV 1 ) X + [SU u Vi] + [U lt 5Vi] = 0, (3.16) 
which can easily be proved. 



Therefore, a binary symmetry constraint of (3.8) can be taken as 

N 

Y a = Z , i.e., [r,^] = [f/ ,$>^ (s ¥ s)T ]- (3-17) 

s=l 

Solving this equation for u, we obtain the required constraints on the potentials 

(3.18) 



U,j = u. 



7* ~ Ij 

where (•, •} denotes the standard inner-product of M. N , B is defined by 

B = diag(/ii,/i 2 ,---,A*iv), (3.19) 

and $i, are defined by 

$ i = (<j> il ,---,<p iN f, ^i = {^n,---,M T , 1<*<4. (3.20) 
The substitution of u with 

U = (&21, "12, "31) "13, "41, "32, "23, "14, "42, "24, "43, U34) 

into the binary Lax systems ( [3.13] ) and (3.14) yields the so-called binary con- 
strained flows: 

4 S > = U(u, X s )^ s \ Vi s) = -U T (u, \ s )iP {s \ l<s<N; (3.21) 

4 S) = V(u,X s )^ s \ 4 S) = -V T (u,X s )^ s \ 1 < s < N; (3.22) 

which are two systems of ordinary differential equations, since V doesn't involve 
any spatial derivative of u. 



3.3 Liouville Integrability 



In order to analyze the Liouville integrability of the above two binary con- 
strained flows, let us first introduce a symplectic structure over 



4 N 



(3.23) 



where £? is defined by (3.19) and <£>,; and \I>j are defined by (3.20). Then the 
associated Poisson bracket reads as 



{f,g} = ^(IdgJdf)=Y / ((g^,B 



i=l 



df D _ x dg v 5/ x dg 



0*i 



4 JV 



-i/ 5/ % 9/ dg 



dipis dcj) ls 84> is dipi 



), f,geC™(R m ), (3.24) 



i=l s=l 

where the vector field /c?/ is uniquely determined by 

uj 2 (X,Idf)=df(X), XeT(R m ). 

A Hamiltonian system with a Hamiltonian H defined over this symplectic 
manifold (R 8JV ,o; 2 ) is given by 

3H QH 
$ it = } = -B" 1 ^, *« = # } = 1 ^ * ^ 4 < ( 3 - 25 ) 

where i is assumed to be the evolution variable. 

In order to prove that two binary constrained flows are Liouville integrable, 
we need generate their Hamiltonian structures and integrals of motion. A 
direct computation can verify the following theorem. 

Theorem 3.1 The spatial binary constrained flow (3.21 ) and the temporal bi- 
nary constrained flow (3.2i) are Hamiltonian systems with the evolution vari- 
ables x, t and the Hamiltonian functions 



k=l 
4 



i<k<i<4 



Oik ~ ui 
Pk -Pi 

0k -0i 
Ik - 1i 



($ fc , (3.26) 



$ fc , £* fc >, (3.27) 



k=l l<fc<Z<4 

respectively, where the matrix A is given by 

A = diag(\i,X 2 , ■ ■ ■ , X N ). 



(3.28) 
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It is well known that Lax representations can engender integrals of motion. 
For the binary constrained flows ( p. 21 ) and ( p. 22 ), let us introduce a Lax 
operator L(X) as follows 

L(X) = (i ii (A)) 4X 4 =T + D(A) =diag( 7l , 72 , 73,74) + D(X), (3.29) 



JV 



1=1 



D(X) = (A,(A)) 4 x4, Ai(A) = ]T^M^, 1 < i,j < 4. (3.30) 
Now if we define 



U(X) = U(u,X), V(X) = V(u,X), 



(3.31) 



then we obtain the following result. 

Theorem 3.2 The binary constrained flows (3.21) and (3.21) have the fol- 
lowing Lax representations: 



(L(X)) X = [U(X),L(X)}, (L(X)) t = [V(X),L(X)}, 



(3.32) 



respectively. 

An r-matrix formulation can also be directly shown for the Lax operator 
L(A) defined by ( jDs| ) and ( |Oo| ). 

Theorem 3.3 The Lax operator L(X) defined by ( 3.2i ) and ( 3.3C ) has the 
following r-matrix formulation 



{L(X)®L(p)} 



1 



fx — X 



V,L 1 (X)+L 2 (n) 



»>i=i 



y Eji, 



(3.33) 



where Li(X) = L(X) <S> I4, -^(m) = I4® L([i), and 

{E V j) M = S ik S jh {L(X) f L(/i)}ii,fci = {Lifc(A),L 3 -j(/i)}, 1 < i,j,k,l< 4. 
Now first from the Lax representations in ( 3.32| ), we obtain 
{vh - i(A)) a - [(7(A), vI A - L{X)}, {vl A - L(X)) t = [V(X), vh - L(X)\, 
where v is a parameter, and thus El 

(det(^/ 4 - L(X))) X = 0, (det(vh - L(X))) t = 0. 
Second from the r-matrix formulation ( |3.33| ), we can obtain 
{tri fe (A),tri'(^)} = 0, k,l>l. 
10 



Expand the determinant of the matrix vl^ — L(X) as 



detail - L(X)) = v A - + F 2 {\)v 2 - Fz{\)v + T 4 (X), (3.34) 

where by Newton's identities on elementary symmetric polynomials, we have 

^i(A) = tri(A), T 2 {\) = |((trL(A)) a - trL 2 (A)), 
.F 3 (A) = i(trL(A)) 3 + ±trL 3 (A) - i(trL(A))trL 2 (A), 
JF 4 (A) = dcti(A) = i(tri(A)) 4 - ±trL 4 (A) + ±(trL 2 (A)) 2 
-i(trL(A)) 2 tr£ 2 (A) + ±(trL(A))tri 3 (A). 

Therefore, it follows that 

{Fi{\)) x = 0, {Ti{\)) t = 0, {^(A),^(/i)} = 0, 1 < i,j < 4. (3.35) 

Make a further expansion of Fi(X) as follows 

= FikX~ k , 1 < i < 4, (3.36) 

k>0 



where obviously the i^o's are all constants. Then it follows from ( 3.3 5| ) that 
(F ik ) x = 0, (Fik) t = 0, {F ik ,Fji} = 0, 1 < i,j < 4, k,l > 0, (3.37) 



Therefore, two binary constrained flows (3.21) and ( 3.2 2j ) have the common 
integrals of motion: Fu, 1 < % < 4, I > 1, which are in involution in pairs 
under the Poisson bracket (3.24). Now it is direct to show the f ollow ing th eorem 
on the Liouville integrability of two binary constrained flows ( |3.21 ) and ( 3.22 ). 



Theorem 3.4 Two binary constrained flows ( 3.21 ) and ( 3.2i ) have the com- 
mon integrals of motion: Fu, 1 < i < 4, I > 1, which are in involution in 
pairs under the Poisson bracket ( 3.24 )> an d °f which the functions Fi S , 1 < 
i < 4, 1 < s < N, are functionally independent over a dense open subset 
of M. SN . Thus, two binary constrained flows fi3~2~l ) and WM) 

are Liouville 

integrable. 



3.4 Involutive Solutions 

Note that u nder the co nstra ints on the potentials J3.18|) , the compatability 
condition of (3.13) and (3.14) is still only the multi-wave interaction equations 
(3.8). Thus, solutions ($i{x, t), ^i(x, t)) to two binary constrained flows ( 3.21 ) 
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and (3.22) will present solutions to the multi-wave interaction equations (3.8), 
namely, 



u ij {x,t) = -± '-(Qifat^BilFjfat)), l<i^i<4, (3.38) 

H ~ Ij 

which also shows the integrability by quadratures for the multi-wave interac- 
tion equations ( |3.8| ) because <&i and ^i, 1 < i < 4, can be determined by 
quadratures. 



It is also easy to observe tha t und er the Poisson bracket (3.24), two Hamil- 
tonians H x and H l defined by ( |3.27 ) and ( 3.27 ) commute, namely 



dH x D _xdH 1 



)-{^B^ dHt 



>) 



(3.39) 



It follows that the above solutions determined by ( 3.38 ) give r ise t o a kind 
of involutive solutions to the multi-wave interaction equations (3.8). Let us 



denote two Hamiltonian flows of ( 3.21 ) and ( 3.22 ) by and gf respectively. 
Then the resulting involutive solutions can be written as 



Uij (^, t) 



7i 

OLi 



■1i 
a 



-{gfg?$> m ,Bgfgfy j0 ) 



J (gf gf$ i0 , Bg* gf* j0 ), 1 < i?j < 4, 



where the initial values <&io and ^io of <I>i and ^i, 1 < i < 4, can be taken to 
be any arbitrary constant vectors of R w . 

In a word, solutions of two binary constrained flows ( [3.21 ) and ( 3^22] ) can 
lead to involutive solutions of the multi-wave interaction equations (|3j). More 
importantly, such involutive solutions show us an explicit Backlund transforma- 
tion from the multi-wave interaction equations ( |3.8| ) to two finite-dimensional 
Liouville integrable Hamiltonian systems determined by ( 3.26 ) and ( |3.27| ) and 
thus the integrability by quadratures for the multi-wave interaction equations 

(H. 



4 Summary and Conclusion 



The multi-wave interaction equations (3 



have been taken as an illustra- 



tive example of binary nonlinearization. Binary symmetry constraints (3.17), 
containing arbitrary distinct constants 71,72,73,74, and arbitrary non-zero 
constants /ii, ^2, M3, /M, were proposed for the multi-wave interaction equa- 
tions (3.8). Two finite-dimensional Liouville integrable Hamiltonian systems 
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determined by ( |3.26| ) and ( p. 27 ), resulted from two binary constrained flows, 



dete rmine a kind of involutive solutions to the multi-wave interaction equations 
(fj|. If we take B = I N , and 

r = V , i.e., diag(7i,72,7 3 ,74) = diag(/3 x , /3 2 , /?3, 

then the result established here will cover all results presented in Ref. E3. 
This also implies that our result leads to a larger class of finite-dimensional 
integrable Hamiltonian systems. 

The theory of binary nonlinearization shows that binary symmetry con- 
straints in the continuous case decompose soliton equations (PDEs) into binary 
constrained flows (ODEs). The r-matrix formulation can be used to expose the 
Liouville integrability of binary constrained flows. The resulting constraints on 
the potentials also present a kind of involutive solutions to soliton equations 
and thus eventually show the integrability by quadratures for soliton equations. 
The whole theory can also be applied to many other types of soliton equations. 

We point out that in the case on 2 x 2 traceless matrix spectral problems, 
binary nonlinearization will be reduced to mono-nonlinearization, if we take 
the reduction 

{i^ij^f = (-0 2J ,<M T , 1 <i < n. 

That is to say, in such a case, all results achieved through mono-nonlinearization 
can be obtained from the results achieved through binary nonlinearization. 
Therefore, in this sense, binary nonlinearization is much broader and more 
systematic. 

We have seen that symmetry constraints yield nonlinear constraints on 
potentials of soliton equations, and put linear spectral problems (linear with 
respect to eigenf unctions) into nonlinear binary constrained flows (nonlinear 
again with respect to eigenfunctions), which seemly makes it more complicated 
to solve soliton equations. However, since spectral problems are overdeter- 
mined, to guarantee the existence of eigenfunctions of spectral problems, one 
needs compatability conditions which are often nonlinear partial differential 
equations. The symmetry property of the constraints brings us the Liouville 
integrability for nonlinear binary constrained flows. Therefore, the symmetry 
property makes up for the complexity of nonlinearization while manipulating 
binary nonlinearization. 

Of particular interest in the study of binary symmetry constraints are 
to create new classical integrable Hamiltonian systems which supplement the 
known classes of classical integrable systems and to expose the integrability by 
quadratures for soliton equations by using binary constrained flows. We men- 
tion that high-order symmetry constraints with involved Lie-Backlund symme- 
tries having non-degenerate Hamiltonians can also be carried out for soliton 
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equations without much difficulty, hut the case of degenerate Hamiltonians 
needs some particular considerations. The rigorous considerations on binary 
symmetry constraints and two-binary symmetry constraints, especially on vari- 
able symplectic structures of binary constrained flows, are being expected to 
be made in a future publication. 
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